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1. Introduction 

The class of TeV-emitting high-mass X-ray binaries (HMXBs) consists of only a few sources: 
LS 5039 [1], PSR B1259-63 [2], LS I +61° 303 [3], HESS J0632-r057 [4], and the newest member 
of the class IFGL J1018.6-5856 [5]. Only the compact object of PSR B1259-63 has been firmly 
identified as a pulsar, while fhe nafures of fhe compacf objecfs in fhe ofher sysfems have nof yef been 
unambiguously defermined. Consequenfly, fhe fundamenfal sefup fhaf produces fhe TeV emission 
along wifh ifs characferisfic variabilify on fhe fimescale of one orbifal period remains unknown. 

LS I -1-61° 303 is fhe only known TeV binary in fhe Norfhem Hemisphere fhaf has a suffi- 
cienfly shorf orbifal period fo allow for regular sfudy wifh TeV insfrumenfs. Locafed af a disfance 
of ~ 2kpc [6], fhe sysfem comprises a BO Ve sfar and a compacf objecf [7, 8]. The observed emis¬ 
sion is variable and modulafed wifh a period of P « 26.5 days, believed fo be associafed wifh fhe 
orbifal sfrucfure of fhe binary sysfem [3, 9-13]. Radial velocify measuremenfs show fhe orbif fo 
be elliptical e = 0.537 ± 0.034, wifh periasfron occurring around phase (j) = 0.275 and apasfron af 
(j) = 0.775 [14]. The periasfron disfance befween fhe sfar and fhe compacf objecf is esfimafed af 
2.84 X 10^° m (0.19 AU) and fhe apasfron disfance af 9.57 x 10'° m (0.64 AU) [15]. The inclination 
of fhe sysfem is expecfed fo lie in fhe range 10° -60° [15], buf fhe lack of exacf knowledge of fhis 
paramefer leads fo considerable uncerfainfy of fhe ofher orbifal paramefers. 

In fhis work, fhe resulfs of fhe VERITAS campaign on ES I -i-61° 303 in late 2014 are presented. 
During fhis time, VERITAS observed hisforically brighf flares from fhe binary around apasfron, 
wifh fhe source exhibiting flux levels a facfor of 2 - 3 fimes higher fhan ever previously observed. 

2. Observations 

VERITAS is an imaging afmospheric-Cherenkov telescope (lACT) array locafed af fhe Fred 
Eawrence Whipple Observatory (FEWO) in soufhern Arizona (1.3 km a.s.L, 31°40’ N, 110°57’ W). 
If consisfs of four 12 m diameter Davies-Coffon design optical telescopes, each wifh a reflector 
composed of 345 fessellafed hexagonal mirror facels fhaf focus lighf onto a camera af fhe focal 
plane wifh a field of view of 3.5°. If is sensitive to photons wifh energies from 85 GeV fo > 30 TeV 
wifh an energy resolution of 15-25%, an angular resolution of < 0.1° af 1 TeV, and fhe abilify 
fo defecf a 1% Crab Nebula source in approximately 25 hours ^ [16]. For a full descripfion of fhe 
hardware componenfs and analysis mefhods utilized by VERITAS, see [10, 17, 18]. 

In fhe 2014 season, VERITAS observations of ES I -i-61° 303 were faken from October 16 
(MJD 56946) fo December 12 (MJD 57003) for a fofal of 24.7 hours of qualify selected livefime. 
These observafions covered fhree separafe orbifal periods, sampling fhe orbifal phase regions of 
0 = 0.5 — 0.2 (see Figure 1 and Table 1). Over fhe entire sef of observafions, a fofal of 449 excess 
evenfs above an energy fhreshold of 300 GeV were defecfed above background, equivalenf fo a 
sfafisfical significance of 21 sfandard deviafions above background (21a, calculafed using Equation 
17 of [19]). 

During fhe firsf orbif observed (in October), fhe source presented fhe largesf of ifs flares (re¬ 
ferred fo as “El” from here on), beginning on 2014 October 17 (MJD 56947, (j) = 0.55) wifh 
emission reaching a peak of (31.9 ± 3.4stat) x 10^^^ cm^^ s^^ on October 18 (MJD 56948). This 
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Orbital phase (j) 


Figure 1: Light curve of LS I +61° 303 during the 2014 observation season in nightly bins. The data 
for the first orbit (October) are shown with orange circles, the data for the second orbit are shown with 
purple diamonds, and the data for the third orbit are shown with blue squares. Flux upper limits at the 99% 
confidence level (using the unbounded approach of [20]) are shown for points with < 3 (J significance and 
are represented by arrows. 


is the largest flux ever detected from this source. Observations following the peak were limited 
by poor weather conditions and only recommenced on October 20 (MJD 56950), by which time 
the flux from the source had already decreased. During the second orbital passage in November, 
VERITAS detected another period of elevated flux (“F2”) from the source at similar orbital phases 
{(j) = 0.5 — 0.6) with peak emission of (20.2 ± 2.2stat) x 10^^^ cm^^ s^' on November 14 (MJD 
56975). 

The rise and fall times of the flares were determined by fitting Equation 7 of [21] to the light 
curve of each orbit and fixing to to the MJD of the observed peak. The fit to the first orbit is not 
very good, with a fit probability of 3 x 10^^. Nevertheless, the rise and fall times of El are found to 
be 0.39 ±0.07 days and 0.37 ±0.23 days, respectively. The second orbit is well fit by this function, 
with a fit probability of 1 x 10^ \ possibly due to the better data sampling throughout this flare. The 
rise and fall times of F2 are thus found to be 0.65 ± 0.13 days and 0.83 ± 0.12 days, respectively. 
A piecewise defined exponential function was also fit to El and E2, but resulted in a very poor fit 
probability for both flares. Variability on a nightly timescale was tested using the method described 
in [22]. Similar to their findings for this source, a hint of nightly variability at a significance level 
of ~ 3 a post trials was found in EL 

Eollow-up observations conducted by VERITAS during the next month (2014 December 10- 
12) covered the orbital phases of 0 = 0.59 — 0.67 and detected the source at a lower flux level, 
reaching only (7.2 ± L2stat) x 10^^^ cm^^ s^^ around the orbital phase at which the flares were 
detected in the previous orbits. The observations during this month appear to exclude the type of 
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Date observed [MJD] 

Orbital phase 

Flux(> 300GeV) [x 10 cm ^ s 

56946 

0.52 

<0.15 

56947 

0.55 

2.32 ± 0.54 

56948 

0.60 

3.20 ± 0.34 

56950 

0.67 

<0.47 

56951 

0.71 

<0.64 

56952 

0.75 

<0.55 

56954 

0.82 

<0.28 

56956 

0.90 

<0.86 

56958 

0.97 

<0.38 

56960 

0.05 

<0.46 

56962 

0.12 

<0.96 

56973 

0.54 

<0.53 

56974 

0.58 

1.30 ±0.34 

56975 

0.61 

2.02 ± 0.22 

56976 

0.65 

1.55 ±0.25 

56977 

0.69 

<0.94 

56979 

0.76 

<0.85 

56981 

0.84 

<0.53 

56982 

0.88 

<1.02 

56983 

0.92 

<0.58 

56985 

0.99 

<0.82 

56987 

0.06 

<0.69 

56989 

0.14 

<1.13 

57001 

0.59 

<0.66 

57002 

0.63 

0.72 ±0.12 

57003 

0.67 

1.06 ±0.20 


Table 1: VERITAS observations of LS I +61° 303 in 2014 


peaked flaring behavior seen at the same phase range in the previous two orbital cycles, perhaps 
indicating some orbit-to-orbit variations in the system. 

The average differential energy spectrum from all observations of LS I +61° 303 during the 
2014 observing season was found to be well fit with a power law of the form 


dN 

dE 



(2.1) 


in which Nq is the normalization at the pivot energy of 1 TeV, and F is the spectral index. The 
measured parameters are consistent with past observations. Differential energy spectra were also 
extracted from FI (October 17-18) and F2 (November 13 -15) and show a similar spectral shape, 
albeit with a higher normalization constant. The parameters from the spectral fits are given in 
Table 2. An uncertainty on the energy scale of 15-25% results in a systematic uncertainty of 
~ 50% on the flux normalization and ~ 40% on the integral flux, assuming a spectral index of -2.5. 


4 







Bright TeV flares from LSI +61° 303 


A. O’Faolain de Bhroithe 


Observations 

Normalization [x 10 cm ^s ^ TeV ^] 

Spectral index 

All (average) 

1.7 ±0.7stat±0.9sys 

—2.35 ± 0.32stat ± 0.3sys 

LI (Oct 17-18) 

8.6±1.0stat±4.3sys 

—2.24 ± 0.12stat ± 0.3sys 

L2 (Nov 13-15) 

4.8 ±0.4stat ±2.4sys 

—2.36 ± 0.12stat ± 0.3sys 


Table 2: Spectral parameters of the power law fits to the observations of LS 1+61° 303 in the energy range 
0.3-20 TeV. 

The systematic uncertainty on the spectral index is estimated at ~ 0.3, accounting for uncertainties 
on the collection efficiency, sky brightness, analysis cuts and simulation model. All spectra are 
shown in Figure 2 along with previous spectral measurements for comparison. 

After the detection of F2 , VERITAS released an ATel [23] notifying the astronomical commu¬ 
nity of the historic flux levels and triggering observations by multiwavelength partners, as well as 
additional observations with the MAGIC TeV experiment. The data from this campaign are under 
analysis and will be presented in an upcoming publication. Initial results of TeV, GeV and X-ray 
observations during this campaign are presented in [24]. 

3. Discussion and Conclusion 

The nature of the compact object in LS I +61° 303 is not firmly established, so proposed 
models of the system cover a range of possibilities. In general, the models fall into one of two main 
categories: a microquasar (ftQ) scenario or a pulsar binary (PB) scenario. The majority of both 
model types employ a leptonic origin of the observed non-thermal emission. 

In the /iQ scenario, non-thermal particle acceleration processes occur in the jet of an accreting 
compact object. Soft photons from the optical star and its accretion disk, and also from the syn¬ 
chrotron radiation in the jet can be upscattered by the ultrarelativistic electrons in the jet [26, 27]. 
Support for this type of model comes from the detection of an elongation of the source on the scale 
of tens of milliarcseconds in radio wavelengths [28, 29], interpreted as a Doppler-boosted jet with 
precession [13, 30]. 

The PB scenario utilizes the presence of a shocked wind in which particle acceleration is the 
result of the interaction between the stellar and the pulsar winds. For example, the author of [31] 
proposes the compact object to be a young pulsar with a spin-down rate of the order of 10^^ erg 
s^* that generates a pulsar wind composed of principally monoenergetic electrons/positrons. The 
pulsar wind expands and shocks as it interacts with the stellar wind from the optical star. This 
shock has a “bow” or “comet” shape with a tail extending away from the optical star. The stellar 
photons can be upscattered to TeV energies by the particles accelerated at the shock. Support for 
this type of model comes from the detection of cometary-shaped radio emission, pointed away 
from the high-mass star [32]. 

A general PB scenario is presented by [33] describing an inhomogeneous stellar wind in which 
the Be star disc is disrupted and fragments, and the resulting clumps of the disc fall into the shock 
region, pushing it closer to the pulsar. The reduction in size of the pulsar wind termination shock 
could allow for increased acceleration efficiency on the timescale of a few hours, depending on 
the size and density of the disc fragments. Such a scenario could qualitatively account for the 
exceptionally bright TeV flares and orbit-to-orbit variations seen in LS I +61° 303. 
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Energy [TeV] 


Figure 2: Average and flare differential energy spectra of LS 1+61° 303 from the VERITAS 2014 observa¬ 
tions, shown in comparison with the average spectra from [10] and [25]. 


The detection of pulsed emission from LS I +61° 303 would unambiguously identify the source 
as a pulsar binary. While no pulsations have been detected to date, it is also possible that the dense 
stellar environment of the source could hinder such a detection. Regardless, further observations 
of LS 1 +61° 303 with TeV instruments are necessary to fully understand the varying TeV emission 
from the source. 
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